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ABSTRACT 

Q\ ■ The next generation of astronomical instruments will be able to detect dense 

pockets of ionized gas created by the first luminous objects. The integrated free-free 
emission from ionized halos creates a spectral distortion of the Cosmic Microwave 
Background, greater than the distortion due to the reionized intergalactic medium. 
Its detection is well within the capabilities of the planned DIMES satellite. Ionized 
halos may also be imaged directly in free-free emission by the Square Kilometer Array. 
Bright halos will be detected as discrete sources, while residual unresolved sources 
can still be detected statistically from temperature fluctuations in maps. Balmer line 
emission from ionized halos is detectable by the Next Generation Space Telescope and 
CnJ ■ can be used to obtain redshift information. Unlike Lya, the Ha line does not suffer 

from resonant scattering by neutral hydrogen, and is still observable when the IGM is 
, not fully reionized. In addition, it is less susceptible to dust extinction. Finally, the 
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kinetic Sunyaev-Zeldovich effect may be used to detect ionized gas at high redshift. 
Used in concert, these observations will probe ionizing emissivity, gas clumping, and 
star formation in the early universe. 



1. Introduction 



^ ' Stimulated in part by the upcoming launches of the next generation of astronomical 

instruments, the epoch of reionization has become the subject of increasingly intense theoretical 
effort in recent years. Reionization affects observations of the Cosmic Microwave Background, both 
as a source of new anisotropics, and in suppressing small-scale fluctuations from the surface of last 
scatter (Haiman & Knox, 1999, and references therein). In addition, the Next Generation Space 
Telescope (NGST) offers the prospect of directly detecting the first luminous objects. Nonetheless, 
the large number of complex, non-linear physical processes involved means that it is extremely 
difficult to make firm theoretical predictions with any kind of confidence. While it appears quite 
likely that NGST will be able to image Population III objects in the infrared continuum (Haiman 
& Loeb 1997, 1998), and detection of the Gunn-Peterson trough may yield the redshift at which 
the IGM is not yet fully reionized (Miralda-Escude & Rees 1998), this information alone falls 
far short of being able to fully constrain models of reionization. Key among the uncertainties is 
the ionizing emissivity of collapsed objects, and the degree of gas clumping. In this paper, we 
suggest the observation of diffuse gas and Population III objects in thermal bremstrahhlung and 
recombination line emission as a direct probe of these quantities. Due to their dependence, 
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free-free and recombination line fluctuations trace the dumpiness and ionization fraction of the 
gas. While the intensities of these signals are directly proportional to one another, their joint 
detection should prove invaluable, given the highly different nature of the foreground contaminants 
and instrumental systematics in each case. In addition, we examine the prospects for detecting 
ionized gas at high redshift via the kinetic Sunyaev-Zeldovich effect. Below, we briefly summarize 
the observational case for each signal. 

Balmer line emission Discussion in the literature on detecting line emission from high 
redshift galaxies has thus far focussed predominantly on Lya (Miralda-Escude & Rees 1998, Kunth 
et al 1998). However, detection in Lya faces problems due to attenuation by dust and resonant 
scattering by neutral hydrogen. Indeed, in the dark ages before the universe was completely 
reionized, it is unlikely that we will be able to detect sources in Lya emission, due to the large 
column density of neutral hydrogen along the line of sight. By contrast, Balmer lines face none 
of these problems, and provide a much cleaner probe of ionizing emissivities. Thus, for instance, 
the intensity of Ha emission is intimately linked to star formation rates. In this paper, we show 
that for equivalent integration times with NGST, moderate resolution spectroscopy (R ~ 1000) 
should yield Ha line detection with the same signal to noise as continuum IR imaging for z < 6.6, 
and 10% signal to noise compared to IR imaging for z > 6.6. Thus, if a high-redshift galaxy is 
detected in imaging mode, Ha line detection should be eminently feasible. In addition, higher 
order Balmer lines such as H/3 should also be detectable. The redshift information provided by 
Balmer line emission should provide an invaluable complement to free-free and IR imaging. 

Free-free emission In addition to Compton scattering off hot gas via the Sunyaev-Zeldovich 
effect, free- free emission from the ionized IGM can also produce a spectral distortion of the 
CMB (Barlett & Stebbins 1991). This distortion increases quadratically at low frequencies, 
and is poorly constrained by COBE observations near the peak of the CMB spectrum {y ~ 56 
GHz). One of the drivers behind the development of the DIMES satellite (Kogut et al 1996; 
see [http:/ /ceylon.gsfc. nasa.gov/DIMES] ) is the hope of detecting this distortion produced by the 
ionized IGM at low frequencies {y ~ 2 GHz). The magnitude of the distortion is directly related 
to the optical depth of ionized gas in the IGM, and thus has the potential to constrain the epoch 
of reionization. However, the ionized IGM is not the only source of free-free emission. Lyman 
limit systems in ionization equilibrium with the external intergalactic ionizing background emit 
thermal bremsstrahlung radiation (Loeb 1996), as do ionized halos with an internal source of 
ionizing radiation (Haiman & Loeb 1997). In this paper, we argue that if the escape fraction 
of ionizing photons from starburst galaxies is small (as is indicated by local observations), then 
the integrated contribution of ionized halos to the free-free background overwhelms that of the 
IGM, regardless of the details of the reionization history of the universe. Thus, the detection of a 
spectral distortion by DIMES is likely to probe the integrated ionizing emissivity and clumping of 
gas in halos, rather than the optical depth of the reionized IGM. 

A natural way to tell if an observed spectral distortion is produced by ionised halos rather 
than an ionised IGM is to observe fluctuations in the free-free background. If the spectral 
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distortion of the CMB is produced by the IGM, it should exhibit a fairly smooth spatial signature 
across the sky, whereas if ionized halos predominate then the free-free background should exhibit 
a great deal of small scale power. Since Galactic foregrounds in free-free and synchrotron emission 
damp rapidly at small scales, it should be possible to directly image the predicted population of 
high redshift ionized halos with a high resolution instrument. This is an application well suited 
for the proposed Square Kilometer Array, or SKA (Braun et al 1998) which will operate in the 
frequency range 0.03 — 20 GHz with an angular resolution comparable to that of HST (~ 0.1") 
and a sensitivity 100 times better than the VLA, down to the ~ 20 nJy level. Such a quantum 
leap in performance should reap immense scientific rewards. 

At present, the appearance of the radio sky at nano Jansky levels is not known. We make a 
first stab at characterising observable properties of a predicted population of free-free emitters. In 
our Press-Schechter based model, we find that ~ 10 4 sources at z > 5 should be directly detectable 
as discrete sources in the 1° SKA field of view. In addition, the integrated background from 
sources too faint to be directly identified creates a fluctuating signal, similar to the Extragalactic 
Background Light (EBL) in the optical and the Cosmic Infrared Background (CIB) in the IR, 
which may be characterised by its power spectrum and skewness. Finally, we investigate the 
clustering of sources at high redshift and find that while it is fairly substantial, the observed signal 
is still likely to be dominated by Poisson fluctuations. 

In all numerical estimates, we assume a background cosmology given by the 'concordance' 
values of Ostriker & Steinhardt (1995): (n m , fi A , Sl b , h, a 8h -i , n) = (0.35,0.65,0.04,0.65,0.87,0.96). 
We use expressions from Carroll, Press & Turner (1992) and Eisenstein (1997) to compute 
cosmological expressions such as the growth factor and luminosity distance. 

2. Modelling the Ionizing Sources 

2.1. Estimating Halo Luminosities 

Since the free-free and Ha emissivity are both proportional to (rig), observational prospects 
hinge crucially on this quantity at high redshift. The signal is dominated by the small-scale 
clumping of the gas at high redshift. One approach to estimating (n 2 e ) is to model the properties of 
halos at high redshift. It is not known to what extent the baryons can cool and condense in halos 
at high redshift. A minimal assumption would be that gas within a halo traces the distribution of 
the dark matter, generally overdense by 8 ~ 180 at collapse; this would indicate (rig) oc (1 + z) 6 . 
Observations of the local universe indicate that the maximum luminosity density of dynamically 
hot stellar systems varies as M -4 / 3 (Burstein et al 1997), which implies that the gas is more 
centrally condensed in low mass systems (which are prevalent at high redshift), presumably 
because dissipational effects are more efficient. Thus, the above scaling with redshift is likely 
an underestimate. If the gas settles to form a disk it could be overdense by 5 ~ 187r 2 A -3 ~ 10 6 
(Dalcanton, Spergel & Summers 1997), where the spin parameter A ~ 0.05 typically (Barnes & 



- 4 - 



Efstathiou 1987). Note, however that (rig) is likely to be dominated by the local clumping of the 
ISM in ionization bounded HII regions (in other words, the clumping factor (n 2 )/(n e ) 2 of the gas 
is significantly greater than unity), in which case these limits would have to be revised upwards. 

A simpler and likely more accurate approach is to assume that most of the ionizing radiation 
produced is absorbed locally in the halo, with only a small fraction escaping (< 20%). This 
is consistent with our observations of the nearby universe. For instance, the observation of 4 
starburst galaxies with the Hopkins Ultraviolet Telescope (HUT) by Leitherer et al (1995) report 
an escape fraction of only 3%, by comparing the observed Lyman continuum flux with that 
predicted from theoretical spectral energy distributions. In our own galaxy, the models of Dove, 
Shull & Ferrara (1999) have an escape fraction for ionizing photons of 6% and 3% (for coeval and 
Gaussian star formation histories respectively) from OB associations in the Milky Way disk, while 
Bland-Hawthorn & Maloney (1999) find an escape fraction of 6% is necessary for consistency with 
the observed line emission from the Magellanic stream and high velocity clouds. The increasing 
density of the gas in halos with collapse redshift indicates that the escape fraction should decrease, 
strengthening our bound. The production rate of recombination line photons -/V" recom b is directly 
related to the production rate of ionizing photons: 

^recomb = 0. B {n\)V « (1 - / esc )iV ion (1) 

where A^; on is the number of ionizing photons emitted per second, and as is the case B 
recombination coefficient. 

In this manner the luminosities of a source in Ha and free-free emission scale directly with the 
production rate of ionizing photons. For Ha, Hummer &: Storey (1987) find that 0.45 Ha photons 
are emitted per Lyman continuum photon over a wide range of nebulosity conditions. We thus 
obtain directly 

L( Ha ) = 1.4 X 10" ergss-. (2) 

For free- free emission, the volume emissivity for a plasma is (Rybicki &: Lightman 1979): 

(rp \ —0.35 
ergs" 1 cm" 3 Hz" 1 (3) 

where we adopt the velocity averaged Gaunt factor gjf = 4.7 and approximate its mild temperature 
dependence with a power law. We can then use the estimate of n 2 V derived from equation (|l]) to 
obtain: 

We thus wish to estimate A^ on (M), the production rate of ionizing photons as a function of 
halo mass. We can either assume that this function is independent of collapse redshift, or that it 
evolves with redshift. In this paper, we adopt the models of Haiman &: Loeb (1998), which assume 
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the former. In particular, we adopt their starburst model, which is normalised to the observed 
metallicity 10~ 3 Z & < Z < 10~ 2 Z & of the IGM at z ~ 3. In this model, a constant fraction of the 
gas mass turns into stars, 1.7% < f s tar — 17% in a starburst which fades after ~ 10 7 years. In 
this paper, we shall assume the upper value Z ~ 10~ 2 Z Q , corresponding to f s t ar ~ 17%, in our 
discussion and plots. We show the results of the calculations for the lower limit star formation 
efficiency f s t ar ~ 1-7% in figure (0) (discussed in section ( [4.2.6D ). During the starburst period, the 
rate of production of ionizing photons is: 

N ion (M) = 2 x 10 53 (/ sto ./0.17) (M/1O 9 M ) photons s" 1 (5) 

Their mini-quasar model yields similar estimates, Ni on (M) = 4 x 10 53 (M/10 9 MQ)photons s _1 , 
with the quasar fading after ~ 10 6 yrs. The quasar model is normalised to the observed quasar 
luminosity function at low redshift, assuming a constant ratio of black hole mass to halo mass. 
We consider the adopted model to be a reasonably conservative estimate for high redshift sources, 
as it seems quite possible that the efficiency of ionizing photon production increases with redshift. 
The star formation rate is tied to the dynamical timescale of the gas, SFRoc Mg/td yn - Since high 
redshift sources are denser, the free-fall time is shorter and iVj on oc SFR oc (1 + Z) 1 ' 5 . Alternatively, 
one can use the empirical Schmidt law SFRoc T, g A (Kennicutt 1998), where T, g is the gas surface 
density, in conjunction with an isothermal sphere model to deduce Ni on oc SFR oc (1 + z) 8 
(all scalings are for an object of fixed collapse mass). Intriguingly, Weedman et al (1998) find 
from observations of galaxies in the HDF that high redshift starbursts have intrinsic UV surface 
brightnesses typically four times brighter than low redshift starbursts. Furthermore, the absence 
of metals to provide cooling in the early universe would tend to lead to a top-heavy IMF (Larson 
1998), further increasing the efficiency of ionizing photon production. Note also that in our 
adopted model the escape fraction should fall with redshift, since A^; on is independent of redshift, 
whereas iV recom b oc (n 2 )V oc (1 + z) 3 (using r v i r oc (1 + z)~ x ). 

We pause to note the possible caveat that supernovae could drive gas out of the galaxy (Dekel 
& Silk 1986), thereby drastically increasing the escape fraction of ionizing photons and vitiating 
our assumptions. The expulsion of gas would also shut off further star formation. Although there 
is enough energy from supernovae to drive such a wind, recent numerical simulations (Mac Low 
& Ferrara 1999) indicate that mass ejection is in fact only efficient in galaxies with M g < 1O 6 M0 
(such galaxies are excluded from our analysis), where the small depth of the potential well means 
that gas is extremely easily unbound. In higher mass galaxies, the supernovae produce a hole 
in the interstellar medium through which the shock heated gas is expelled, while the cold gas 
remains bound. We thus ignore the effect of gas expulsion by supernovae. Note also that our 
emissivities are calibrated to a quiescent, steady-state mode of star formation, as is observed in 
the disk of many spiral galaxies today. There also exists the possibility of massive starbursts, with 
star formation rates of up to 100 M & yr^ 1 , in which a large fraction of the gas mass is converted 
into stars in the short dynamical time of the dense galactic center. This is particularly likely to 
happen as halos merge, and would produce significantly brighter objects, easily detectable with 
the proposed observations. 
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When computing surface brightnesses, it is necessary to assume a characteristic size for an 
ionized halo r. We can assume that the gas profile extends out to the virial radius (r ~ r v i r ) 
or that the gas is confined to a disk with spin parameter A (r ~ Ar„j r ), where the virial radius 
r vir = v c /(10y/2H(z)) = 2.1Kpc(M/lO 9 M ) 1 / 3 (l + z/10) -1 , for an isothermal sphere. In 
subsequent calculations we shall assume r ~ Xr v i r . The small corresponding angular size means 
that in general we do not resolve the galaxy, which may thus be treated as a point source. We 
quote angular size as 6 = max ( telescope , r/dA(z)), where ^telescope is the angular resolution of 
the telescope. Note that this assumption only affects our surface brightness calculations, since 
Lua, free-free oc N ion , independent of r. 



2.2. Statistics of Halos 

We employ a Press-Schechter based model of the abundance of ionizing halos similar to 
that adopted by Haiman & Loeb (1998). In our simplified physical model, halos collapse, form 
a starburst lasting t Q = 10 7 years, and then recombine and no longer contribute to the free-free 
background (FFB). We ignore sources ionized by an external background radiation field, since the 
limiting column density before these sources become self shielding corresponds to a undetectably 
small flux at high redshift. We adopt a different expression for the halo collapse from Haiman 
& Loeb (1998); their expression slightly underestimates the halo formation rate as it includes a 
negative contribution from merging halos. The more accurate expression we use is given by Sasaki 
(1994): 

dNf™ ,^ ldDdnps 51 

-^AT^ Z) = 15~dT~dM~ a^JdjLfi (6) 

where D{z) is the growth factor, and <5 C = 1.7 is the threshold above which mass fluctuations 
collapse (note that at high redshift, the difference between this expression and the one used 

Iff on 
dW 



by Haiman & Loeb (1998) is negligible, since mergers are unimportant). Here, dN I M — is the 



collapse rate of halos per mass interval per unit comoving volume, while ^jjf- is the standard 
Press-Schechter number density of collapsed halos per mass interval (Press & Schechter 1974, 
Bond et al 1991). Then, given an expression for the expected flux from a halo of mass M at 
redshift z, S = S(M,z), we can calculate the expected comoving number density of ionized halos 
in a given flux interval as a function of redshift: 

dN halo rt(z)-t dN^ orm dM 

-dSdV^ z)= J t{z) dt ^tM~l[S (?) 

where t Q is the duration of a starburst. 

Note, however, that star formation or quasar activity is only able to proceed if gas is 
able to cool and condense in the halo. Thus, we set a cut-off mass for a halo to be ionized of 
M* = 10 8 (1 + z/10)~ 3 / 2 Mq, which is the critical mass needed to attain a virial temperature of 10 4 
K to excite atomic hydrogen cooling. With this in mind, we can compute moments of the intensity 
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distribution due to sources above a redshift z m j n via: 

(S n (> z mm ,S c )) = r dz f SmaX dS^-^-S n (8) 
Jz mm Js mi M) dSdV dzdn 

where S c is a flux limit above which discrete source identification and removal is possible. We shall 
use this expression entensively in ensuing sections. The zeroth moment of equation (S) can be used 
to compute the number counts of sources above the flux limit S c . In this case, S max — > oo and 
Smin(z) = max(£ c , S*(z)), where we denote S*(z) as the flux from a halo of minimum mass M* 
at redshift z. One can also use equation (||) to compute the moments of the residual background 
signal due to sources below the flux limit S c , when all discrete sources above the flux limit have 
been identified and removed. In this case, we set S max = S c and S m i n {z) = S*(z). 



2.3. Comparing Halo Emissivity with the IGM 

Throughout this paper we only compute the contribution of ionized halos to free-free and 
recombination line emission, ignoring the emission from filaments and the IGM. While it is obvious 
that small scale fluctuations will be dominated by ionized halos, it is not immediately obvious 
that the mean sky-averaged signal (as is probed by spectral distortion measurements of the CMB) 
from ionized halos is greater than that from other sources. This in fact is a direct consequence of 
our assumption that the escape fraction of ionising photons from halos is small. This implies that 
the degradation of ionizing photons into recombination line photons take place largely in halos. 
The fraction of ionizing photons produced which are degraded in halos as opposed to the IGM 
(including ionized filaments and Lyman limit systems) scales as (1 — (f e sc))/{fesc), where f esc is 
averaged over all halos. This ratio is large due to the small assumed / esc , which (as we've argued 
above) also is likely to decrease with redshift. We can compare the volume averaged emissivities of 
the IGM and ionized halos more quantitatively. The comoving free-free emissivity of the ionized 
halos is given by: 

/•oo f l AT 

^) = l^M—(M.,)L A M,,) (9) 

where L u (M, z) is given by equations (|j) and (||) . The comoving free- free emissivity of the 
IGM may be estimated directly from equation (]3|), assuming the entire IGM is ionized. For 
the purpose of this calculation, we have assumed a uniform IGM and ignored the contribution 
of overdense structures without an internal source of ionizing radiation. This is a good 
approximation as the low ambient radiation field in the early universe means that these objects 
become self-shielding at low column densities and are the last objects to become reionized 
(Miralda-Escude, Haehnelt, & Rees, 1998). In Figure (|l]) we display the relative halo-IGM 
emissivity e = (e) halo /( e ) IGM = (^e)haio/(^e)iGM- For z < 20, the halo emissivity dominates over 
the IGM emissivity (at high redshift, we underestimate e as the IGM is only partially reionized. 
As computed by Haiman & Loeb (1998), full reionization in this model only occurs at z ~ 13). 
Note that this ratio also applies directly to recombination line emission. We conclude that the 
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free-free and recombination line emission is dominated by discrete sources rather than a diffuse 
background. 



3. Recombination line emission 
3.1. Ha emission 

The Ha flux density from an ionized halo as detected with a R=1000 filter may be estimated 

as: 

J - = l^i^K^HiM)(i^) nJy (10) 

where we have estimated Ln a from equations (|2|) and (|5|). Note that the line width of Ha is 
simply given by Doppler broadening AA/A ~ 10 _4 (6/30kms _1 ) (unlike the Lya line, which will 
have broad damping wings), so it is still not resolved at R=1000. Using a narrow bandpass 
enables us to reduce the amount of background noise. Let us compare with the prospects for 
observability in the infra-red, which corresponds to rest-frame UV emission. The UV continuum 
emission in the rest frame 1500-2800 A(longward of the Lyman limit) for a Salpeter IMF is 
L v = 8 x 10 27 (n^ppr) = 8 x 10 27 ( <W ^^^t) ergs^Hz" 1 (Madau 1998). Using equation 
(|5|), this translates into an observed IR flux density: 

Lv , x fl + z\- 1 ( M \ 



Andl y ' \ 10 J V1O 9 M 0/ 

Hence, Jhq/^r ~ 40(12/1000). Note that the spectral energy distribution/IMF of a luminous 
source can be probed by its color, while the ratio of the Ha luminosity to the rest frame UV 
luminosity will constrain the escape fraction of ionizing photons. 

Let us examine the competitiveness of Ha with continuum IR imaging with NGST. The signal 
to noise ratio of an observation with exposure time t is given by (Gillett & Mountain 1998): 

S Id 



(12) 



N (I s t + I bgt + n-I dc -t + n- N?) 1 / 2 

where I s = 1.5 • 10 7 vl • e • 1/R ■ j v ■ / electrons s _1 is the signal photocurrent (A is the area 
of the primary mirror in m 2 , e is the detector quantum efficiency (e > 80% for A < 5/xm, 
e > 50% for A > 5/xm), j u is the source strength in Jy, and f is the fraction of source photons 
collected (/ « 0.7 for A < 3.5^m and / « 0.5 for A > 3.5/iin)), I bg = 1.5 • 10 7 A ■ e ■ 1/R ■ A, • Q 
electrons s -1 is the background photocurrent (/?„ is the background in Jy sr _1 ), n ~ 4 is 
the number of pixels required to cover Q O.Olarcsec 2 , Idc is the detector dark current, 
and N r < 15eread _1 (Stockman et al 1997) is the read noise per pixel. We assume a point 
source subtends Q w max(l, (A/3. 5/xm)) 0.01 arcsec 2 . For A = 0.5 — 5/im, the projected dark 
current is Id c < 0.02es _1 for a Si:As IBC array, while for A > 5/um with an InSb array, 
Idc < les -1 (Stockman et al 1997). To estimate the sky background, we consider twice the 
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sky background observed by COBE (Hauser et al 1998). The sky background is measured at 
1.25, 2.2, 3.5, 4.9, 12, 25, 60/iin. We assume the full wavelength dependence of the sky background 
in our calculations, employing a spline interpolation between the observed points. Roughly 
speaking, the background is < 10~ 5 Jy arcsec -2 in the range A ~ 1 — 5/rni, and rises to 
~ 5 x 10" 4 Jy arcsec -2 for A > 10 fim. An additional source of noise is thermal emission from 
the primary mirror, which we approximate as a 75K blackbody (Barsony, 1999). This exceeds the 
sky background at A ~ 10/iin, and thus only affects sources in Ha emission with z > 15. For IR 
continuum imaging (R = 3) at A ~ 1 — 5/im, the sky background is the dominant source of noise. 

Thus, we have S/N(IR) ps 0.5(^/0.01 arcsec 2 )" 1 / 2 (t$t)~ (lo^fe) fl/2 - So > for instance, a lOcr 
detection of our canonical 10 9 M© source at z=9 would take 400s. For spectroscopic applications 
with R > 1000, the detector noise is the dominant source of noise. The increase in dark current 
at A > 5/mi corresponds to reduced a signal to noise ratio for sources with z > 6.6, leading to 
S/N(Ha;z < 6.6, z > 6.6) « (0.7, 0.06) (17/0.01 arcsec 2 )" 1 / 2 (t^) 1 (losfe) tV2 - A 10(7 detection 
of our canonical 10 9 M Q source at z=9 in Ha would take 3 x 10 4 s. Note that since the flux density 
oc R, the signal to noise is independent of the spectral resolution R in the dark current limited 
regime R > 1000. Thus, even higher resolution spectroscopy, up to the point when the Ha line is 
resolved (R ~ 10 ), should be possible. We derive the useful rule of thumb that: 

<6.6„> 6.6) = <n, 0.!2) (£) (13, 

Thus, for sources which sit at z <6.6, a source may be detected in IR imaging and Ha spectroscopy 
in equal amounts of time. For sources at z>6.6, when the detector dark current for Ha detection 
is signficantly larger, detection in Ha requires an integration time 100 times longer than in IR 
imaging to achieve the same signal to noise ratio. This ratio will be reduced if advances in 
technology permit smaller dark currents at A > 5/im. 

Note that the detection rate depends on the angular extent of the sources, even in the dark 
current limited regime. An extended source covers more pixels, resulting in a higher dark current 
(in very extended sources, it may be best to focus on the central pixels, which contain most of the 
light). For this reason, we quote the relative signal to noise in equation (|l~3|), which is independent 
of the angular size. In figure (Q), we show the number of sources detectable with a S/N of 10 
with R > 1000 spectroscopy in 10 4 s (~ 3 hours) of integration time, all computed using the full 
expression (p^), for both disk-like (r ~ Ar) and extended (r ~ r„j r ) objects. At least several 
thousand objects with z > 5 may potentially be detected in the NGST 4' x 4' field of view. The 
actual detection rate, of course, depends on the targeting strategy for spectroscopy (see discussion 
below). Note that if the star formation efficiency is lower, then the detection rate falls. In figure 
©, we show the detection rate for f star ~ 1.7%, rather than f star ~ 17%. For a star formation 
efficiency lower by an order of magnitude, an integration time 100 times longer is required to 
obtain the same number counts. 



It should be possible to observe higher order Balmer lines as well. For case B recombination 
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with n e = 100cm 3 (the density dependence is fairly weak), /h/?//ho = 0.35 (Osterbrock 1989). 
Accounting also for the different bv = v/R for Ha and H/3, we find that S/N(Rf3)/S/NQla) = 0.25 
for a given source, which implies an integration time 16 times longer is required to observe a source 
in H/3 with the same signal to noise ratio as in Ha (note that we assume that attenuation effects 
are still unimportant for H/3, which is of shorter wavelength than Ha and thus more suspectible 
to dust extinction or stellar absorption). However, in the redshift range 6.6 < z < 9.5, the H/3 line 
falls in the regime A < 5//m with low dark current, while the Ha line falls in the A > 5[im regime 
with high dark current. In this case S ' j 'N (H/3) / S j 'N (Ha) = 2, i.e., the H/3 line is potentially easier 
to observe than the Ha line. Observing the H/3 line should provide a firmer redshift identification, 
while the Ha/H/3 ratio gives a useful measure of extinction. In figure (0), we show the detection 
rate of the H/3 line as a function of redshift. 

How can we pick out potential high redshift targets for spectroscopy? They may be quickly 
identified via continuum imaging using the Lyman-break technique successfully used at low 
redshifts (Steidel et al 1996). Beyond the redshift of reionization, an analogous technique using the 
Gunn-Peterson trough to detect dropouts can be used. Photometric redshifts may be supplemented 
by a narrow band filter search (R ~ 100 — 400) to isolate Ha bright galaxies in a narrow redshift 
range. The filter width involves some trade-offs: the filter should be broad enough to cover a 
sufficiently large comoving volume, yet narrow enough that the flux density in Ha greatly exceeds 
continuum flux densities. Note that current plans for NGST allow for multi-object spectroscopy 
in the NIR 1 — 5/jm range, using micro-mirrors but only single slit spectroscopy in the MIR range, 
5 — 20/zm (Stockman et al 1997). However, even in the NIR a grating spectrograph will be able to 
perform spectroscopy for only a few hundred (~ several percent) of the objects in the field of view 
at a given time, thus requiring many pointings to complete a survey. Alternatively, a promising 
proposal is the Infrared Fourier Transform Spectrograph (IFTS) (Graham et al 1998) which will 
be able to perform panchromatic observations in the 1 — 15/xm range (which is the desired z < 20 
range for Ha observations). It will be able to acquire both broad band imaging data and higher 
spectral resolution data (R=l-10000) for all objects in the field of view simultaneously, with very 
high (near 100%) throughput. Thus, the need for object preselection is obviated, and such an 
instrument would be extremely well suited for our purposes. 

3.2. Comparison with Lya 

While ionizing objects are expected to be more luminous in Lya than in Ha, detection 
in Ha has a number of nice features. Firstly, since Ha is emitted in rest-frame optical, it is 
less susceptible to attenuation by dust than Lya, or rest-frame UV continuum. It is thus a 
significantly cleaner tracer of the star formation rate. Furthermore, the resonant scattering of 
Lya photons within the HII region and in surrounding HI regions means that they will diffuse 
over a large region, resulting in an unobservably low surface brightness. The large optical path 
of Lya photons as they random walk out of these regions greatly increases their probability of 
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absorption by dust (Chariot & Fall 1993). Finally, the observability of Lya is affected by complex 
velocity flows of the neutral gas, which may be outflowing from the ionizing source (thus, for 
instance, the observation of P-Cygni type profiles in local star-forming galaxies (see Kunth el 
al 1998 and references therein)). Note that 50% of the objects selected by the Lyman-break 
technique show no emission in Lya, while the rest have weak rest- frame equivalent widths < 20 A 
(Steidel et al 1996). By comparison, the expected rest-frame equivalent widths of star-forming 
galaxies from modelling their spectral energy distribution is W Q ~ 100 - 200 A (Chariot & Fall 
1993). Also, note that quasars are expected to have equivalent widths in the quasar rest frame 
of W a = 0.68^ j™dvf u /u « 827a~ 1 (3/4) a A = 360A(a = 1.5) (Chariot & Fall 1993), where 
vl, v a are the Lyman limit and Lya frequencies, f c is the continuum flux near Lya, and f u oc v~ a 
blueward of Lya. In fact, most bright quasars are observed to have rest-frame equivalent widths 
of 50 < W a < 150A (Schneider, Schmidt & Gunn 1991, Baldwin, Wampler k Gaskell 1989). This 
could be due to covering of the AGN source by HI clouds or dust extinction effects. We tentatively 
conclude that even in the post-reionization epoch, a Lya survey is an inefficient means to search 
for ionizing sources. Until recently, Ha surveys of high-z objects have not been carried out due to 
the difficulties of infrared spectroscopy from the ground. It is intriguing to note that recent such 
Ha surveys (Glazebrook et al 1998) show a star formation rate approximately three times as high 
as that inferred from the 2800 A continuum luminosity by Madau et al (1996). This reinforces 
the expectation that surveys in UV continuum and Lya tend to underestimate the star formation 
rates as they are subject to highly uncertain dust extinction effects. 

Most importantly for our purposes, a Ha survey would be able to detect ionizing sources 
prior to the epoch when HII regions overlap, when the universe was optically thick to Lya 
photons. In order for Lya emission from a source to be unaffected by resonant scattering, the 
ionizing source must ionize a region large enough for the Lya photon to redshift out of the 
Gunn-Peterson damping wing. This problem has been considered by Miralda-Escude (1998). 
Given that the optical depth of the Gunn-Peterson damping wing as a function of displacement 
AA from line center, r(AA) = 1.7 x 10~ 3 (1 + z/10) 3 / 2 (AA/A) _1 , a region r ~ 1 Mpc in proper 
size must be ionized, independent of redshift, in order to have t(AA) ~ 0.3. Note that this 
is somewhat greater than the proper radius of influence of a single ionizing source at z=9, 
Influence ~ 0.1Mpc(M/10 9 M Q ) 1 / 3 (l + z/10) -2 (assuming the source never reaches its equilibrium 
Stromgren radius). Nonetheless, the clustering of ionizing sources means that many should exist 
in large, overlapping ionized regions. Even if an ionized region r ~ 1 Mpc is created, the high 
recombination rate due to the higher IGM densities are high redshift results in a residual neutral 
fraction which is optically thick to Lya. If we take the column density to be Njji ~ xn^r, where 
x is the ionization fraction in ionization equilibrium, then the optical depth to Lya scattering in 
transversing this region is: 

where J21 is the background ionizing intensity in units of 10 _21 ergs _1 cm~ 2 Hz -1 sr _1 , and 
$b = 8pb/pb is the gas overdensity. It therefore seems dubious that unscattered Lya can be 
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observed except in highly underdense regions, or unless ionizing sources are considerably more 
luminous than we have considered. 

Loeb & Rybicki (1999) have considered the observational signature of Lya emission of an 
ionizing halo, given that the Lya photons are likely to resonantly scatter in the surrounding 
neutral IGM, before redshifting out of resonance. They find that the Lya photons emitted by a 
typical source at z ~ 10 scatter over a characteristic angular radius of ~ 15" around the source 
and compose a line which is broadened and redshifted by ~ 10 3 kms _1 . However, the large angular 
scale leads to a very low surface brightness. Adopting the same estimates as the previous section 
for instrumental noise and the IR background, we find that the signal is likely to be overwhelmed 
by instrumental noise. Assuming A^ya ~ Aj on , so that the flux /L ya ~ 12/hq and the flux density 
jhya ~ 2.2jH Q (Case B), the relative signal to noise ratio at the same spectral resolution is given 

(in fact, line broadening for Lya is likely to further lower this ratio by a factor of a few). Thus, 

with S/N(Lya) = 6 x KT^ft/TOOarcsec 2 )- 1 / 2 (to 2 ) ^ (io»fe) * 1/2 ' our canonical source 1O 9 M 
source at z=9 would require an unacceptably long integration time of 32 days for a 10<r detection. 
If even trace amounts of dust are present, the large number of resonant scatterings means that 
the Lya signal is likely to be severely attenuated. We conclude that this would be a very difficult 
observation, although extremely interesting if carried out successfully. In particular, it would be 
very interesting to view the same bright source in Lya and Ha. Their relative intensities would 
establish the importance of dust extinction, while their relative spatial extent would constrain 
the topology of neutral gas surrounding the source. Note that Loeb & Rybicki (1999) assume a 
uniform IGM surrounding the source. If the surrounding IGM is inhomogeneous, the Lya photons 
are likely to diffuse along the underdense directions, leading to a signal of even lower surface 
brightness. In addition, Lyman limit systems will absorb a certain fraction of the Lya photons 
(since photons scattering in these systems are no longer redshifting out of resonance). 



4. Free-free emission 

Free-free emission from galaxies forms a fluctuating background in radio frequencies similar to 
the EBL sought in the optical (Vogeley 1997) and the CIB sought after in the infra-red (Kashlinsky 
et al 1996, Dwek et al 1998). It is has the potential to yield a great deal of useful information 
in that it is a direct tracer of star formation and clumping, it is unaffected by dust, there are no 
unknown K-corrections to be applied (due to the flat spectrum nature of free- free emission), and 
sources at very high redshifts (z > 10) make a significant contribution and are still potentially 
observable. We adopt two approaches in characterizing the Free- Free Background (FFB): detecting 
the mean signal (via spectral distortion of the CMB) and detecting the fluctuating signal (by 
direct imaging of ionized halos or detecting temperature fluctuations across the sky) . 
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4.1. Spectral distortion 

Consider first the mean sky averaged signal (S). The integrated emission from ionized halos 
creates a spectral distortion of the CMB which is potentially observable. Note that for a given 
mean surface brightness (S), the temperature perturbation rises quadratically at low frequencies, 
where (in the Rayleigh- Jeans limit) ATff = c 2 (S)/2kv 2 . This spectral distortion may be 
characterised by the optical depth to free-free emission (Bartlett & Stebbins 1991) Yff = T ff x 2 
where T 7 is the undistorted CMB photon temperature, x = hv/kT-y and AT// = T e fj — T 7 is the 
mean temperature perturbation created by free-free emission (note that the Y// has a different 
spectral form from Compton y-distortion (parameterized by y), arising from Compton scattering 
of CMB photons off hot gas). The quadratic nature of the distortion means that a bound on 
Yff is best obtained by comparing temperature measurements at high frequencies (when the 
perturbation due to free-free emission is small) and low frequencies (when it becomes important). 
The CMB blackbody temperature near the peak of emission has been measured by the FIRAS 
instrument aboard COBE to be T 7 = 2.768 ± 0.004K (95% CL) (Fixsen et al 1996), while at low 
frequencies the effective temperature has been measured as T e g = 2.651q;|qK at 1.4 GHz (Staggs 
et al 1996a) and T eff = 2.730 ± 0.014K at 10.7 GHz (Staggs et al 1996b). Combining the limits on 
extant data to compute AT//, one obtains Y ff < 1.9 x 10~ 5 (95% CL) (Smoot & Scott 1996). 
Significant amounts of free-free emission can still be present without violating this constraint. At 
low frequencies, it corresponds to a limit on the mean temperature perturbation: 

44 * io ~ 2 k < i6 > 

One of the chief goals of the DIMES satellite (Kogut 1996) is to measure spectral distortions 
from free-free emission at low frequencies. DIMES will have a frequency coverage 2-100 GHz and 
a sensitivity O.lmK (at 2 GHz). It is hoped that a measurement of AT// will indicate the optical 
depth of the ionized IGM and thus constrain the redshift of reionization z r . However, we have 
already seen from quite general arguments in section 2.2 that the mean emission from ionized 
halos should dominate that of the IGM. One can examine this in more detail. The equation of 
cosmological radiative transfer (Peebles 1993) gives the observed surface brightness of the ionized 
IGM as: 

1 f z r dl e(v(l + z)) 
Air Jo dz (1 + z) 

where e is given by equation (||). If z r = 13, as in our model, then this translates into 
AT// = c 2 J/2ku 2 = 6 x 10" 6 K at 2GHz, which is beyond the O.lmK sensitivity of DIMES, which 
will only be able to detect z r > 85 for the ACDM cosmology we have adopted. On the other 
hand, if we compute (S) from equation (Q) (setting z m - m , S c = since point source removal is not 
possible with the 10° beam of DIMES), we obtain AT// = 3.4 x 10~ 3 K, well within the capability 
of DIMES, while still satisfying the y-distortion constraint by an order of magnitude. Note that 
we have assumed that free-free emitters are the only radio-bright sources in the sky. The existence 



j(u ) = ^i ^ t: ,T,3 (17) 



-14- 



of radio-loud galaxies and AGNs(see section 4.3.3 for discussion) is likely to increase the spectral 
distortion due to radio sources. 

Thus, it is likely that any spectral distortion detected by DIMES will not constrain the 
reionization of the IGM but rather the integrated free-free and other radio emission from galaxies, 
which provides a stronger signal by at least an order of magnitude. 



4.2. Direct Detection of Ionized Halos 

A natural way to determine if an observed spectral distortion is due to free-free emission from 
the first ionized halos or the IGM is to attempt to directly detect sources with a high angular 
resolution instrument. 



4-2.1. Detection of point sources 
The free- free flux from an ionized halo may be estimated as: 

5„ = « (1 + z) „ 2.5 (l±i\ " ( " ) (*) '°- 35 „ Jy (18) 

4vr4 v ' v io J vio 9 M y vio 4 k; j v ; 

where we have estimated from equations (||) and (|5|). Assuming a point source, and the 0.1" 
resolution of the Square Kilometer Array (SKA), this translates into a brightness temperature of: 

* - — (^r w (wk) m l (^r <-> 

On the other hand, the detector noise may be estimated as (Rohlfs & Wilson 1996): 

2kT SYS ( Au \ - 1 / 2 ( t \ - 1 ' 2 , , 

where we have used ^4 e fj/T sys = 2 x 10 8 cm 2 /K for the SKA (Braun et al 1998), and we assume a 
bandwidth Av w 0.5u. Thus, in 10 days, a 5a detection of a 16 nJy source is possible. This is a 
tremendous leap in sensitivity; the deepest observations with the VLA at 8.4 GHz (Partridge et 
al, 1997) with approximately a week's total integration time were able to identify point sources at 
the 7 fiJy level (the rms sensitivity was 1.5/iJy), with a resolution of 6". 

For a given integration time, the sensitivity for detection of a given halo in free- free emission 
with the SKA is significantly less than that for Ha detection (with R > 1000) with NGST. The 
relative signal to noise is: 

f (free -free), x , , 9x / Au \ 1 ' 2 f t s \ 1 / 2 

i'( ! <7,o7) = (l.lxl0- 1 .3xl»-)(— ) (£) (21) 
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where the above ratio holds in the regime where free-free detection is limited by instrumental 
noise (for a discussion of confusion noise, see section ( 4.2.5D ). Thus, the SKA will only be able 
to detect bright sources. However, note that its field of view is considerably larger than that of 
NGST-1 rather than 4'. In figure (|3|) we compute the number of sources detectable above a given 
flux threshold S c - We find that ~ 10 sources with z > 5 should be present in the 1 square degree 
field of view of the Square Kilometer Array above a source detection threshold of 70nJy. 



4-2.2. Power spectrum of unresolved sources 

Identification of individual sources may not always be feasible. Many sources will too faint 
to be individually detected, or they may be crowded and blended together in a confused field. 
Such sources can still be detected statistically by observing the fluctuations in the FFB due to 
sources below the flux limit. The discreteness of sources create large amounts of small scale power. 
The power spectrum is white noise (C; ~ constant) until the angular scale of the beam (or the 
typical scale of a halo, if it is larger), beyond which it damps rapidly. We can compute compute 
the statistics of the background below the flux cut S c , assuming the sources are unclustered. In 
particular, we compute the power spectrum of the fluctuations from equation (|8]): 

c Poisson = ^ = / " S 2 dS (22) 

Jo oo 



(independent of I), and thus compute the rms temperature fluctuations from 
T rms = J2i ~4^rCi lsson Bf « \- K e' 1 where we have taken the beam of the telescope to be 
Gaussian: Bf = e~ l ( l+1 ^ 9 . Note that the rms temperature fluctuations increase with increasing 
resolution. In Figure (Q) we display the computed rms temperature fluctuations at 2 GHz and 
with 4" resolution. Also shown is the mean signal T, computed from the first moment of equation 
(^|), which is independent of angular resolution. Our model is consistent with present observations. 
With S c = 7^Jy and = 6" we obtain AT/T = 9.5 x 1(T 5 at 8.4GHz, as compared with the 
measurement AT/T ~ 7 ± 8 x 10 -5 of Partridge et al (1997), and within their 95% confidence 
limit AT/T < 1.3 x 10" 4 . 



Asssuming that galactic foregrounds are negligible (see section ( 4.2.5[ )), the main source of 



noise is instrumental. From equation (|20|), the rms sensitivity to temperature fluctuations for a 
beam of angular size 9 is: 

T - = « lt U ) M (iGHi) [WS) (23> 

For the same parameters, and assuming a cut-off flux of S c = 70nJy, our expected signal is 
only 0.02K. However, note that our expected signal T rms oc 0~ l , whereas the instrumental noise 
2~instrum OC 6 . By varying the size of the interferometer, we can vary the angular size of the beam, 
optimising it for this observation. Referring to figure (||), we see that on angular scales ~ 1" — 10", 
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the temperature fluctuations due to the free-free background exceeds both the instrumental 
noise level and the galactic synchrotron foreground. Therefore, observing a radio field at coarser 
resolution increases the signal to noise ratio of temperature fluctuations sufficiently to enable the 
fluctuating free-free background to be detected. Note that if we simply smooth the map, the 
?lnstrum oc 6~ 2 scaling is not strictly applicable. Smoothing the map corresponds to downweighting 
longer baseline information, which essentially means that some information is thrown away. This 
decreases the effective area and increases the detector noise, so Ti nstrum declines more gradually 
with 6 (although still faster than T rms oc 9~ l ). The best solution would be to observe in a more 
compact configuration with resolution ~ 1" — 10" to begin with, so no data is thrown away. 
One wants to reach a compromise between the need for the high angular resolution required 
to beat down source confusion in discrete source identification and removal, and the somewhat 
lower angular resolution optimal for analysing the fluctuating residual signal, or confusion noise 
itself. A natural solution would be to work at an angular resolution # e quai where instrumental and 
confusion noise are comparable. Working at higher resolutions than this does not yield significant 
advantages for discrete source extraction, since the instrumental noise background in flux units 
is independent of angular resolution. By fiat, after source removal the instrumental noise and 
fluctuating signal are now about equal, and smoothing the fluctuating signal increases the signal 
to noise ratio sufficiently for an analysis to be undertaken. Note that source removal need not 
always be in the instrument noise dominated regime, S c — n& ^instrum* Source identification in a 
crowded field can be difficult, and not all sources can be removed down to the flux limit. If the 
angular size of typical objects is large, 6>t yp > e quai> then confusion noise exceeds instrumental 
noise and S c > n CT S'i ns trum- In this case, increasing the angular resolution beyond 9 typ yields no 
advantage. Both because of the higher cut-off flux S c , as well as the lower angular resolution, 
detecting the fluctuating signal will not be a problem. Furthermore, as we discuss in section 
( 4.2.5| ), our model probably underestimates the temperature fluctuations which will be observed, 
which is likely to have components from other radio souces in the sky such as low-redshift radio 
galaxies and radio-loud AGNs. 

The free-free component may be separated from other contaminants in multi-frequency 
observations by its spectral signature (flat spectrum, as opposed (for instance) to the B v oc i/~ ' 9 
spectrum of synchrotron radiation). In addition, besides merely computing RMS temperature 
fluctuations, it would be desirable to compute the angular power spectrum of the observed 
signal and confirm its white noise character. In performing this analysis one has to take care to 
ensure that there is no aliasing of power on scales larger than the field of view (due to Galactic 
foregrounds and the CMB). This can be done by high-pass filtering and edge tapering the map. 
Besides the various sources of noise depicted in figure ([]), one also has to contend with sample 
variance, given by: 



where / s k y is the fraction of sky covered (this expression assumes the C\ 's are Gaussian distributed 
as they should be by the Central Limit Theorem, since a large number of modes contribute). In 




(24) 
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addition, one probably will bin the computed C/'s in bins of size Al to increase the signal to 
noise. Using 18 ~ n, an angular resolution of 9 ~ 4" corresponds to / ~ 10 5 . The signal to noise 

for the sample variance term is -^r = 20 (j^) (l^Ti) (wo) ' wn ^ cn ma y ^ e increased by- 
increasing the bin size Al. Thus, sample variance should not pose any problems. 



4-2.3. Clustering of ionized sources 

Does clustering of the ionized halos make a significant contribution to fluctuations in the 
FFB? We find that on arcsecond scales and with a cutoff flux S c = 70nJy, Poisson fluctuations are 
more important than the contribution due to clustering. However, with a lower cutoff flux S c , or 
on somewhat larger angular scales, the clustering term dominates. 

First, let us consider the evolution of the halo correlation function with redshift. The matter 
correlation function may be computed as: 

Um{r, z) = r ^4irk 3 P(k, z)™^ (25) 
jo K kr 

However, halos which collapse early in the history of the universe consitute rare density peaks 
which are highly biased. In a linear bias model, the correlation function of halos is given by: 

£ hh (M 1 ,M 2 ,r,z) = b(M 1 ,z)b(M 2 ,z)Um(r,z) (26) 

where the bias factor is given by Mo & White (1996): 

b(M, z) = l+ V —-± = 1 + - 1 (27) 

where 5 C = 1.68 and v = 5 c /a(M, z). For simplicity we assume the correlation function of halos 
has the same shape at different redshifts, and differs only in amplitude. This assumption is shown 
to hold in numerical simulations (Kravtsov & Klypin 1998, Ma 1999), due to a cancellation of 
non-linear effects in the evolution of the power spectrum and bias factor. We assume £ mm (r) is 
a power law on scales below the non-linear lengthscale (since this is shown to hold in numerical 
simulations), and is given by equation ( pq ) on scales above the non- linear lengthscale. In Figure 
(|5|), we compute the halo correlation length r^ al ° (defined as Chh(^o alo ,z) = 1) at each epoch z, 
assuming a number weighted bias factor 

SS.KM,z)n(M,z)dM 

» {z) ~ TSMKm (28> 

The cutoff mass M* = 10 8 (1 + z/10)~ 3 / 2 M@ introduces strong bias at high redshift, with the net 
result that the correlation length of halos decreases very slowly with redshift, despite the fall in 
the matter correlation length. For a higher cutoff mass, bias becomes important at lower redshift, 
and the initial dip seen in Figure (||) moves to lower redshift and is less pronounced. Thus, the 
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observed clustering in the Lyman-break galaxy sample of Giavalisco et al (1998), with a comoving 
correlation length at z ~ 3 comparable to that in the local universe, is consistent with a detection 
of the most luminous and massive (M ~ 8 x 10 11 h~ 1 MQ) galaxies at that redshift. The clustering 
of ionising sources at high redshift is therefore non- negligible, and a potentially important effect 
in contributing to free- free background fluctuations. 

The angular correlation function of the free-free background below a flux cut-off S c is: 

0( , lft) = f'yr^SMr^,,, (29 ) 
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where j e ff(S c ,z) is the proper free-free emissivity of halos at redshift z below the flux limit 
S c , and b(S,z) = b(M(S,z),z). We use the small angle approximation r 2 = A 2 + d\9 2 , an 
excellent approximation since the clustering length covers a small part of the sky, r Q /dA <C 1. The 
contribution to the angular power spectrum due to clustering cf nster is then straightforwardly 
computed from the Legendre transform of the correlation function: 

P /"OO 

^duster = 2qi J d(cos9)P l {cos9)C{6) « 2n J 9d9J {l9)C{9) (30) 

where the last step is valid for I 3> 1. 

In Figure (0), we show the dependence of the clustering term on the flux cutoff at fixed 
angular scale 9 ~ 1". As we lower the flux cutoff S c , the Poisson term declines more sharply than 
the clustering term, until the clustering term dominates. We can understand this intuitively: 
the Poisson term is dominated by rare bright objects just below the detection threshold. As we 
lower the detection threshold, such objects are excluded from the residual FFB, leaving a more 
uniform distribution of faint objects. By contrast, the clustering term, like the mean brightness 
term, is dominated by faint, low mass objects. In particular, the amplitude of the clustering term 
is determined by the bias of objects at the mass cut-off, M c = 10 8 (1 + z/10)~ 3/2 M Q , rather than 
the bias of the rarer (though more highly biased) high mass objects. Thus, subtraction of bright 
sources decreases the Poisson term much more effectively than the clustering term. In Figure (0), 
we compute the amplitude of the clustering contribution to temperature fluctuations as a function 
of angular scale at fixed flux cutoff S c = 70nJy. We find that on the arcsecond scales on which the 
surface brightness of free-free fluctations exceeds instrumental noise and Galactic foregrounds (see 
below), the Poisson term is more important than the clustering term. However, at a lower flux 
cutoff the clustering term would dominate on these scales. Note the different angular dependence 
of the two terms: Cf msson ~const, whereas clustering ^ w h ere £ hh ^ r ~7. On large scales, 

the matter correlation function falls below a power law, and eventually Cf uster cx l n (where 
the linear power spectrum P(k) oc k n ). This downturn is unimportant on the scales plotted, 
though it becomes appreciable on larger scales. At such low fluctuation amplitudes, both galactic 
foregrounds and intrinsic CMB anistropies dominate over the free-free signal. 
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4-2.4- Skewness 

Another signature of residual discrete sources is the skewness they induce in maps 
(Refregier, Spergel & Herbig 1998). We compute this as the third moment of equation 
(§), h = ^ p ixjQ C dS^S 3 Jy 3 . The instrumental noise is Gaussian distributed, and thus 
has vanishing skewness. The la measurement uncertainty on this vanishing skewness is 
v[Vinstrum\ = Jj^vfnstrumi where (Jinstrum, is the rms noise per pixel. Thus, the signal to 

noise of a skewness measurement is jj = f^fp) ^ ( s Is — )• We assume (Jinstrum ~ 1/^Jy 

V / Y instrum J 

for the VLA, and (Jinstrum ~ lOnJy for the SKA, and also that point sources may be detected 
and removed at S c = 7/iJy for the VLA and S c = 70nJy for the SKA. The number of 
pixels is given by N pix = (FOV/9 pix ) 2 , where FOV=40' for the VLA and FOV=l° for the 
SKA. We assume pixel sizes of order ~ 10" for the VLA and ~ 1" for the SKA. From 
equation (|8|) we obtain:^ = (1//Jy) 3 for the VLA and ^3 = (12nJy) 3 per pixel for the 
SKA. This translates into S/N ps 5 3 (6 /W")(a- mstvum /lfiJy) 3 (S c /7 fiJy) 1 - 4 -^ for the VLA and 
S/N « 13 3 (0/l")(a instruin /lOnJy) 3 (5 c /7OnJy)( 4 -^ for the SKA, where |§ cx S' 13 in the 
neighbourhood of S c . The skewness is thus a very efficient way of detecting the residual field, 
potentially much more powerful than attempting to use the second moment. It uses the fact that 
free- free sources can only contribute positive temperature fluctuations, whereas fluctuations due 
to instrumental noise are symmetric about zero. It is thus a fairly robust statistic. Depending 
on how accurately Gaussian the instrumental noise is, the same exercise can be carried out with 
kurtosis. Note that we have assumed that Galactic foregrounds have no skewness. Given that 
their second moment damps rapidly at small scales (see below), both the amplitude and variance 
of the third moment is likely to be small. It would be interesting to look at existing VLA maps 
for residual skewness. 



4-2.5. Foreground contamination 

A possible concern is foreground source contamination. Note that due to the finite thickness 
of the surface of last scattering, CMB fluctuations damp rapidly beyond I ~ W 3 (9 ~ 0.1°). 
Thus, intrinsic CMB fluctations from the surface of last scattering contribute negligibly on the 
scales probed here, Cf MB ~ 0. One can filter out the large scale power with a filter function 
F\ = 1 — e^ +l ^ a f ilter , where cr/ater ~ 0.1°. At small scales, the Galactic foregrounds C\ oc B{y)l~ 3 
should damp rapidly(Tegmark & Efstathiou 1996, Bersanelli et al 1996), although note that the 
empirical l~ 3 power law fit has only been observed for I < 300 (Wright 1998); we assume this 
behaviour may be extrapolated to smaller scales. Below, we briefly summarise the estimated 
contributions from various foregrounds. 



Galactic Free-free emission This component is dangerous as its spectral signature is 
identical to our signal. However, the large scale component of the Galactic free-free emission can 
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be removed by cross-correlation with Ha maps. For instance, the WhaM survey (Reynolds et al 
1995; see |http:// www . astro . wise . edu / wham| ) will map the northern sky at 1° resolution down to 



0.01 Rayleighs (1 Rayleigh=2.41 x lO^ergs^cm^sr- 1 ), or T ff = UfiKT^ - 35 (^Hi) 2 ' 15 - We 
assume B{v) oc j/ -015 and normalise to the DIRBE analysis of Kogut et al (1996), which obtained 
rms fluctuations of 7.1 [iK at 53 GHz (COBE's finite beam of course causes an underestimate of 
the rms fluctuations, but because C\ oc Z -3 this is a negligible effect). This yields: 



c free-}ree = ^ ^ 1Q _ 4 /~ 3 R 2 S r (31) 

which translates into 5Tf ree ~f ree = 2.3 x 10~ 6 {^hiY 2 '^ K - 

Galactic Synchrotron We assume B{y) oc i/~ 0,9 and normalise to an rms fluctuation of 11 
\iK at 31 GHz, the upper limit of Kogut et al (1996), as indicated by a cross-correlation between 
408 MHz and 19 GHz emission (de Oliveira-Costa et al 1998). This yields: 

Cr = 6.1 x 10- 3 (^) T 3 K 2 sr (32) 

which yields ST syn = 8.6 x 10~ 6 (qT 7 ) ( 2GHz ) 2 9 ^ ne synchrotron emission may be identified 
and eliminated by its distinct spectral signature. 

Extragalactic Point Sources Extragalactic point sources are essentially the signal we are 
looking for. In this paper, we predict the number of sources in a field of view above a certain 
flux S c which may be directly detected, as well as attempt statistical characterizations of the 
fluctuating background signal due to faint undetected sources. Here we pause to note two facts. 
Firstly, there are many other sources in the radio sky (radio galaxies, AGNs, etc) than the 
free-free emitters we have attempted to characterize. Secondly, when detecting point sources, the 
fluctuating signal due to undetected sources is also a source of confusion noise. 

How do the figures in our model compare with existing observations? The deep radio survey 
by Partridge et al (1997) obtain a maximum likelihood estimate of the integral source count as: 

/ a \ -1.2±0.2 

N(> S) = (17 ± 2) j arcmin- 2 (33) 

Vl/uJy/ 

Assuming that most sources have flat spectra l v oc v a , where a ~ in this frequency range, we 
see from figure (|3|) that our model yields fewer sources than the observed number counts at 7/iJy 
by a factor of a few. Our admittedly crude model is not meant to be accurate at low redshift and 
ignores many contributions to the radio sky. Our low source counts, together with the fact that we 
satisfy present y-distortion bounds by an order of magnitude, indicate that the model employed 
is very conservative and probably underestimates the observable signal. On the other hand, 
extrapolating the observed source counts to fainter flux levels clearly overestimates the signal, as 
the mean flux (S) diverges and violates the y-distortion constraint. Evidently the source counts 
must turn over at low flux levels (in our model, this occurs because of the low mass cutoff M*). 
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To be conservative, we underestimate the signal to noise ratio of projected observations. Thus, 
we employ our semi-analytic model to estimate the fluctuating background signal, whereas when 
estimating confusion noise in point source detectation, we use the extrapolated source counts. 

Note that we are not completely powerless to distinguish against contaminants, such as low 
redshift AGNs. If star formation occurs in an extended fashion across galactic halos, we should be 
able to spatially resolve our target sources, whereas AGNs would only show up as point sources. 
In addition, by observing at several frequencies, the observed spectral index a (where S v oc u~ a ) 
can be used to diagnose the origin of radio emission. In particular, it can be used to distinguish 
between thermal and non-thermal emission. Finally, as a zeroth order cut, our target high redshift 
galaxies should be amongst the faintest objects in a sample, both because of increased cosmological 
dimming, and the fact that collapsed objects at high redshift are intrinsically less massive and 
luminous. 

We now estimate the confusion noise. Using equation (|33]), we can place an upper bound on 
the noise from residual undetected sources from (S 2 ) = J^ c j^S 2 dS = 25.5 5°' 8 //Jy 2 arcmin~ 2 , 
which yields the white noise power spectrum: 

/ O \ 0.8 / x -4+2a 

This gives 5T PS = 0.09K 1 [j^)° A (^dk) ^' b y far the leading contribution 

to foreground contamination at these frequencies and angular scales. We can compare the 
relative contribution of instrumental and confusion noise by demanding that the source 
flux cutoff be n a ~ 7 times the rms noise. In this case, using equation (p0[), we obtain 

fe^f = L (b40 (¥)°- 8 {mt) ' 5 (ukf 5 (M*- In the confusion noise dominated regime, 
the appropriate cut-off flux may be determined by directly comparing the flux received from a 
source at the cut-off flux S c and the background flux below this cut-off: 

s_ „ s c = f^_\ - 6 (JL\ 1 

N ~ ((5 2 }| Sc ^)V2 UonJyy 1 {o.4» J { > 

Confusion noise thus becomes the leading source of noise in discrete source identification and 
extraction at low angular resolution or if sources are extended rather than point-like, i.e., 
9 = maX^bcamAbjcct) > 0.4". 

Will source crowding be a problem? If we extrapolate the Partridge et al (1997) counts down 
to 70 nJy, then there will be ~ 400 sources per square arcmin, or ~ 9arcsec 2 per source. With a 
maximum resolution of 0.1", source crowding of bright point sources such as disks or mini-quasars 
should not be a problem. Even if the sources are extended, the angular extent of typical high 
redshift M* objects is small, 9 v { r ~ r v i r /dA < 1". However, it is important to note that nothing is 
known about the appearance of the radio sky below 1/xJy. It is possible that a host of extended 
low redshift, low surface brightness objects may suddenly surface at these flux levels. 
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4-2.6. Lower Star Formation Efficiency Case 

Let us consider the lower star formation efficiency case where only 1.7%, rather than 17%, 
of the gas mass in halos fragments to form stars. This corresponds to an IGM metallicity at 
z=3 of Z ~ 10~ 3 Zq, rather than Z ~ 10~ 2 Z & . The net effect of the lower value is that while 
the number of sources remain the same, each source is an order of magnitude fainter. Thus, 
to obtain the same number counts as previous plots, an integration time 100 times longer is 
required. Also, the brightness of the fluctuating background is reduced, and is equivalent to 
removing discrete sources an order of magnitude fainter. The net result is that point source 
detection is likely to be instrument noise limited rather than confusion noise limited. Thus, 
after discrete source removal, the fluctuating residual background will be significantly harder 
to detect. In figure (0), we show results for the lower star formation efficiency, assuming 
the same integration time and source removal threshold as before. Note from the plot of 
foregrounds that the fluctuating background signal is only marginally detectable. Besides 
the RMS signal, the skewness is also reduced, to ^3 = (0.5/iJy) 3 for the VLA (down from 
^3 = (l^Jy) 3 ) and ^3 = (5.5nJy) 3 for the SKA (down from ^3 = (12nJy) 3 ). Skewness at 
this level is still detectable, with S/N=2.5 3 (9/W")(a instTum /lfiJy) 3 (S c /7 fi3y)^-^ (VLA) and 
S/N=5 3 (0/l")(a instrum /lOnJy) 3 (5 c /7OnJy)( 4 -^ (SKA). 

The situation in reality is likely lie between the two scenarios we have considered, and thus 
the actual number counts is probably bracketed by our calculations. The temperature fluctuations 
in the residual field are detectable if discrete source removal is confusion limited; otherwise, it 
seems likely that only the non-gaussian signature of the residual field will yield useful information. 

5. Kinetic Sunyaev-Zeldovich effect 

Sunyaev-Zeldovich effects arise from the scattering of CMB photons off moving electrons, 
which result in brightness fluctuations of the CMB. They are a very attractive means of detecting 
high redshift objects, in part because of the well-known fact that the surface brightness of 
SZ effects is independent of redshift. This is because the increased energy density (1 + z) 4 of 
the CMB at high redshift keeps pace with surface brightness dimming effects. Thus, since the 
angular size of an object of fixed physical size increases with redshift beyond z ~ 2, the SZ 
flux from a given object increases with redshift. The kinetic effect is given by Tcmb = c r 
while the thermal effect is given by (fcT/m e c 2 )r, where the optical depth r = J n e aTdr. For 
the photoionized halos we are considering, typical temperatures are T ~ 10 4 K. Their peculiar 
velocities D(z) 1 ' 2 kms 1 , where D(z) is the growth factor (normalised to 1 today) 

and the characteristic one-dimensional velocity dispersion of galaxies in the local universe 
is v ~ GOOkms- 1 (Strauss & Willick 1995). Thus, (v/c)/(kT/m e c 2 ) ~ 400(1 + z/10)" 1 / 2 
and the kinetic effect dominates over the thermal effect. For an isothermal halo, we obtain 

1O 9m ) K> similar to the temperature fluctuation produced by a 
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cluster in kinetic SZ at low redshift (the high redshift halo is ~ 1000 times smaller, but ~ 1000 
times denser). 

Unfortunately, this is likely to be difficult to detect. The flux from a halo is related to the 
temperature shift by: 

2k 3 T 2 r 2fc 3 T 3 v N 

S u = -^q(x) j dn\AT u (e)\ = _^Lg( x )-crr^ (36) 

where q(x) = x 4 /4sinh 2 (x/2) is the spectral function, ,x = hu/kT^, T 7 = 2.7K is the CMB 
temperature, and N e = J n e d 3 r ~ M/m p (Qb/^m) is the total number of free electrons in the 
system. From the form of the spectral function, the kinetic SZ effect peaks at v = 217 GHz. This 
is well beyond the frequency coverage of the SKA, which has an upper limit of 20 GHz. At this 
frequency, and taking v ~ v rms , we obtain a flux density S u (20GHz) = 0.24( 10 ^ fQ )(l+-z/10) L5 nJy. 



Comparing with equation (|i8|), the SZ flux density at 20 GHz is ~ 10 times less than the free-free 
flux density. We can maximise the SZ flux by going to higher frequencies. At the peak frequency, 
S u (2l7GHz) = 9.6( 1Q9 A f / )(1 + z/10) L5 nJy, and the kinetic SZ signal dominates over free-free 



emission. At these frequencies, the MMA (see http: / / www.mma.nrao.edu ) would be the premier 



instrument, with a resolution of 0.1" at 230 GHz. However, its rms sensitivity at this frequency is 
0.21 mJky/min ' 5 , which means that even in a week's integration time it will only be able to go 
down to //Jy levels. We conclude that SZ detection of high-redshift ionized halos is not possible in 
the near future. 

However, it is worth noting that the SKA and MMA will be able to detect the kinetic 
SZ effect from the HII regions in the IGM blown by the first luminous sources. If the escape 
fraction is f esc ~ 20% for our 1O 9 M source, it will be able to create an ionized region with 
N e ~ Ni on t ~ 2 x 10 67 electrons (using equation (|5|) and assuming the source lasts for t Q ~ 10 7 yrs). 
This corresponds to a region r comov i n g = (3A e /47rn como „ in9 ) 1 / 3 ~ 1.2 Mpc in size (note that 
Nion/Vp TopeT aB "-proper ~ 10~ 7 <^ 1, so the effect of recombinations is negligible). Since this is 
much smaller than the coherence length of the velocity field, we assume all parts of the bubble 
are moving with a uniform peculiar velocity. From (|36|) , we find that the flux from such an 
ionized bubble is S' l /(20GHz) = 35nJy (detectable by the SKA, which, using equation ( p0| ) and 
A eff /T sys = 1 x 10 8 cm 2 /K at 20 GHz (Braun et al 1998) has a Si nstrum — 6nJy rms sensitivity) 
and S U {217GH z) = 1.4 /^Jy (detectable by the MMA). However, along any given line of sight 
one intersects not one but a multitude of ionized bubbles. Thus, the signal has to be detected 
statistically. Calculating the CMB fluctuations due to inhomogeneous reionization has been 
the subject of much recent work (Agahanim et al 1996, Grusinov & Hu 1998, Knox et al 1998; 
see also Haiman & Knox 1999, Natarajan &; Sigurdsson 1999), and is beyond the scope of this 
paper. However, we make two observations. Firstly, predictions have focussed on the possibility 
of detection by Planck (detection by MAP seems unlikely). The temperature fluctuation power 
spectrum is generally white noise with a peak on the scale of a typical bubble; the induced 
AT/T ~ 10~ 6 — 10~ 7 . For the anistropies to be detectable by Planck, the bubbles need to 
have an unrealistically large size r com oving ~ lOMpc. For bubbles with r comov i ng ~ IMpc, the 
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anisotropies peak on sub-arcminute scales, beyond the angular resolution of Planck (6>fwhm ~ 5.5' 
at 217 GHz; see bttp: / / astro.estec.esa.nl/Planck/ ) . High resolution ground based instruments 
such as the SKA or MMA, which can easily resolve the bubbles, might be better suited to the 
task. Furthermore, the significantly larger collecting area of these instruments means that they 
will be more sensitive than Planck by about 3 orders of magnitude: with a year of observing 
time, the rms sensitivity of Planck at 217 GHz is at best ~ lmJy. If we require that the sample 
variance be sufficiently low that 5 = C//AC/ = 10, then using equation (|2~I|), we find that a map 
of size = 2(5/10) 2 (//10 4 ) _1 (A//10 2 ) -1 square degrees is sufficient. This is easily achievable 
with the SKA. Secondly, the signal due to inhomogeneous reionization (IR) is quite possibly 
smaller than that due to the Ostriker-Vishniac effect, which arises from the coupling between 
density and velocity fluctuations in a fully reionized medium (Haiman & Knox 1999, Jaffe & 
Kamionkowski 1998). To isolate the IR signal, it would be interesting to cross-correlate observed 
CMB fluctuations with the observed distribution of ionizing halos. They could either correlate 
positively (if halos blow ionized bubbles around themselves) or negatively (if underdense voids are 
reionized first). 



6. Conclusions 

In this paper, we have used a simple model based on Press-Schechter theory to make 
observational predictions for the detection of high redshift galaxies in free-free and Ha emission. 
These signals are valuable as they directly trace the distribution of dense ionised gas at high 
redshift. Furthermore, they are relatively unaffected by dust or resonant scattering, and their 
intensities are directly related to one another. They can thus serve as a very powerful probe of the 
dark ages, and constitute promising applications of the new generation of instruments designed 
to probe this epoch, in particular the Next Generation Space Telescope (NGST) and the Square 
Kilometer Array (SKA). 

A fairly robust prediction of our model is that the integrated free- free and other radio emission 
from ionizing sources will swamp the free-free emission from the ionized IGM by an least an order 
of magnitude. Thus, any spectral distortion of the CMB observed by DIMES at low frequencies 
will constrain the former quantity. This claim depends only upon the small escape fraction of 
ionizing photons from sources, as is observed in galaxies in the local universe. To evade it, the 
ionizing sources must themselves contain very little ionized gas, either because gas has been driven 
out by supernovae, or ionizing photons are able to escape through a hole in the ISM, leaving most 
of the gas in the host galaxy neutral. 

The escape fraction of ionizing photons f esc from high redshift halos may be constrained by 
comparing the rest frame UV luminosities longward of the Lyman limit of the sources as observed 
by NGST, with the expected luminosities in Ha and free-free emission. Given a model for the 
relative spectral energy distribution, the relative luminosities then directly constrain the escape 
fraction in halos at high redshift, an important quantity in theories on reionization. If this fraction 
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is low, as we argue, then in a given amount of integration time, the Ha emission from an ionized 
halo should be detected (with R=1000) with the same signal to noise as IR continuum imaging 
with NGST for z < 6.6. For z > 6.6, due to the increase in dark current noise in the A > 5fim 
regime, it will be detected with 10% of signal to noise ratio of IR imaging. Higher order lines such 
as H/3 should also be detectable. In addition, the SKA should be able to directly detect the same 
ionizing sources in free-free emission. While the NGST will be able to detect fainter sources than 
the SKA, the SKA has a larger field of view: ~200 times larger in solid angle. We find that the 
SKA should be able to detect ~ 10 4 individual free-free emission sources with z > 5 in its field of 
view. However, a large population of other radio bright sources (AGNs at low redshift, etc) will 
also be detected. NGST (which can obtain redshift information from the Ha line) must therefore 
be used in conjunction with the SKA to distinguish between these two populations. 

Free-free sources below the flux limit of the SKA constitute a fluctuating background 
which can be detected statistically. At arcsecond scales, Poisson fluctuations in the number of 
free-free sources create temperature fluctuations detectable by SKA. At such small scales, galactic 
foregrounds are expected to damp rapidly. We find that the expected clustering of high-redshift 
ionizing sources makes a negligible impact on the observed power spectrum of temperature 
fluctuations. These temperature fluctuations have a highly non-Gaussian signature. In particular, 
they should result in a detectable skewness. Finally, ionized gas at high redshift may also be 
detected by the kinetic Sunyaev-Zeldovich effect. While individual ionized halos are probably too 
faint to be detectable by the SKA or MMA, the ionized IGM around these halos create CMB 
anisotropies on sub-arcminute scales which should be detectable by these instruments. 

I am grateful to my advisor, David Spergel, for his encouragement and guidance. I also 
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Fig. 1. — The relative halo/IGM free-free emissivity e = (e)h a io/(e)iGM as a function of redshift, 
assuming that the entire uniform IGM is ionized. The contribution of halos to the comoving free- free 
emissivity dominates that of the IGM at all redshifts below 20. 
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Fig. 2. — Number of sources in 4' x 4' NGST field of field with a redshift greater than z, detectable 
with Ha, H/3 spectroscopy (solid and dashed lines respectively) as a 10a fluctuation with t = 10 4 s, 
and R=1000, for emission from a disc (r ~ Xr v i r , essentially a point source), and an extended source 
(r ~ Tyi^*) • The dip at z — 7,9 for Ha, H^/3 is due to increased detector noise at longer wavelengths 
A > 5 /im. Note therefore the existence of an interval where H/? is potentially easier to detect than 
Ha. 
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Fig. 3. — Number of sources which may be detected in the 1° field of view of the Square Kilometer 
Array (the field of view for the VLA is similiar, 40' x 40'), as a function of the cut-off flux S c . 
Realistic limiting fluxes for point source detection are shown. The extrapolated source counts from 
Partridge et al (1997) are also shown. 
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Fig. 4. — Mean and RMS temperature fluctuations for diffuse background below point source 
removal threshold S c , at 2 GHz and 4" resolution. Note the steeper dependence of Poisson 
fluctuations on the flux cutoff S c compared to the clustering term. 
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Fig. 5. — Characteristic number weighted correlation length r Q (comoving) as function of redshift, 
and corresponding angular scale O = rP ropcr /c?A- The bias of high redshift halos means that angular 
clustering does not plummet with redshift. For a higher cutoff mass M#, bias becomes important 
at lower redshift. Beyond the turnoff, the angular scale of clustering is almost constant at a few 
arcminutes. 
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Fig. 6. — Angular dependence of T rms = \ l<yl+ ^ Cl } 1 ^ 2 at 2 GHz, as a function of the beam FWHM 
9. Here the FFB is composed of sources with S < S c = 70nJy. Temperature fluctuations of the 
FFB increase with higher angular resolution, whereas the foregrounds damp rapidly. Also shown 
is the instrumental noise for integration times of 1 and 10 days. There exists a window in the 
2" — 20" range where the signal is greater than both instrumental noise and foregrounds. Note that 
the Poisson term is greater than the clustering term in this domain. 



-35- 



Ha number counts 



~i — i — i — i — i — i — i — i — i — i — | — i — i — r 

Ha 

H|8 




J I I 1_J I I I L_l IXJ I L 



6 8 10 12 

Z 

Free-free background 
-2.5 



E- 1 

O 



i i i i i i i i i i i i i i i i i i i I- 



<T2>l/2 

Poisson 




< T 2 > l/2 " 
^ A cluster: 



i i i i i i i i i i i i i i i i i i i i- 



-8 -7 -6 -5 
log S„ (Jy) 



Free-free number counts 
6 




-8 -7 -6 -5 
log S c (Jy) 
Free — free foregrounds 

-2 



W "3 



Z -4 



o -5 



-6 - 




0.1 1 10 100 1000 
#(arcsec) 



Fig. 7. — Figures (Q) (Ha number counts), (||) (Free-free number counts), (01) (Free-free 
background), (||) (Free- free foregrounds), computed for the lower star formation efficiency case 
fstar = 1-7%. All other assumptions are unchanged. 



